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We investigate the procedures of loading and detecting three-dimensional fermionic quantum gases 
in a one-dimensional optical superlattice potential subjected to a trapping potential. Additionally, 
we consider the relaxation dynamics after a sudden change of the superlattice potential. We numer¬ 
ically simulate the time-dependent evolution of the continuous system using exact diagonalization 
of non-interacting fermions. During the loading procedure we analyze the occupation of the instan¬ 
taneous energy levels and compare the situation in a homogeneous system with the trapped one. 

Strong differences are found in particular in the evolution of excitations which we trace back to the 
distinct global density distribution. Starting from an imbalanced state in the superlattice potential, 
we consider the relaxation dynamics of fermions after a slow change of the superlattice potential 
and find a bimodule distribution of excitations. To be able to compare with the experimental results 
we also simulate the measurement sequence of the even and odd local density and find a strong de¬ 
pendence of the outcome on the actual ramp procedure. We suggest how the loading and detecting 
procedure can be optimized. 


I. INTRODUCTION 

The recent experimental progress to manipulate quan¬ 
tum gases made it possible to study different interesting 
quantum phenomena. The clarification of the mecha¬ 
nism underlying the BCS to Bose-Einstein condensate 
crossover in different dimensions has been one of the suc¬ 
cessful advances p}. Further, in optical lattices the su- 
pcrfluid to Mott-insulator transition in bosonic gases [2} 
and the liquid to Mott-insulator transition in fermionic 
gases 0 H} have been realized. However, a lot of in¬ 
teresting quantum phases such as the antiferromagnet, 
d-wave superconductivity or spin liquids M have not 
yet been reached due to the low temperature or more pre¬ 
cisely low entropy per particle required [f| 3]. For bosons 
it has been shown that the temperature increases with 
ramping up an optical lattice and it can be minimized 
with tuning the trap frequency §, f|. Thus, the detailed 
investigation of the loading procedure of a fermionic gas 
in three-dimension is of utmost importance. 

Moreover, quantum gases have the advantage that they 
are only very weakly coupled to the environment and 
can be considered as almost isolated on relatively long 
time-scales. Thus, quantum gases enable the investiga¬ 
tion of the quantum dynamics of isolated systems. A 
high amount of attention have received so-called quan¬ 
tum quenches. A quantum quench is one specific way to 
set the system out of equilibrium: starting with a sys¬ 
tem prepared in an equilibrium state of a chosen initial 
Hamiltonian, one or several of the system parameters are 
changed. The state evolves according to the new Hamil¬ 
tonian. Such a change is either sudden or with a certain 
ramp speed. The quench velocity of the slow quench 
plays an important role and has been studied in gapped 
and gapless phases and across phase transitions relating 
to the Kibble Zurek mechanism 0 - One of the fun¬ 
damental questions related to a sudden quench is how 
fast correlations can be formed and whether the isolated 


system relaxes to a thermal equilibrium state. These 
questions have been widely studied theoretically both in 
integrable and non-integrable systems as in Refs. flQH38| . 

So far many experiments were conducted for bosonic 
gases. The observation of the dynamics of bosons 
in optical lattices [39l - l46| , the prethermalization of 
one-dimensional condensates [43, 0 , the light-cone-like 
spreading of correlations after a quench [is} are examples 
of the experimental activities in this area. Moreover, the 
dynamics within a superlattice has been studied in the 
experiment [50| and theory 

Less experiments have been conducted for fermionic 
quantum gases. Mainly, interaction or lattice quenches 
JM1] and spin dynamics have been investigated [57|, 

In this paper we study the ramp up procedure to 
load fermionic atoms into first two-dimensional pancakes 
and then into a one-dimensional superlattice potential. 
The additional remaining two spatial dimensions are sub¬ 
jected to a harmonic trap. We analyse in detail the load¬ 
ing procedure and suggest improvements for the exper¬ 
iments. Starting from a prepared imbalanced quantum 
gas, we perform a fast and slow quantum quench and fol¬ 
low the relaxation dynamics of the fermionic gas in the 
trap. In particular, we find a bimodule density distribu¬ 
tions at long times which can be attributed to a situa¬ 
tion in which locally the gas is equilibrated, but globally 
a highly excited state is reached. Additionally, we anal¬ 
yse in detail the detection sequence of the local even-odd 
density and point out possible improvements to the se¬ 
quence. 

II. DESCRIPTION OF THE SYSTEM 

We study a non-interacting gas of N = 10 5 (for ex¬ 
ample 40 K) fermionic atoms with mass m in a three di¬ 
mensional dipole trap and a one-dimensional optical su¬ 
perlattice along the ^-direction. The Hamiltonian of the 
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non-interacting atoms at time t is given by 
H= I dr 


WVr)— —V 2 ty(r) 
w 2 m w 


+ \rn{ul{t)x 2 + u 2 y {t)y 2 + u 2 {t)z 2 )^ 

+ F(x,f)'I' t (r)'I'(r) . (1) 


Here 'l't(r) and 'F(r) are the fermionic creation and an¬ 
nihilation operators. The bichromatic lattice potential 
along x-direction is formed by an infra-red laser beam 
with wavelength A R = and a green laser beam with 

wavelength Ac = At time t the potential can 

be described by the superposition of the two lattices 

V(x,t) = Vc(t) cos 2 (2k R x - - V R (t) cos 2 (fc_Rx)(2) 


where V R {t) and Veit) are the amplitudes of the red and 
green lattice potentials at time t , respectively. The ge¬ 
ometry of the optical superlattice is strongly determined 
by the phase difference (f> between the two lattices at the 
position of the atomic cloud. In the experiment the phase 
can be controlled by detuning the frequency of the green 
laser from exactly twice the frequency of the red laser. 

At the position of the atoms the dipole trap can be 
approximated by a harmonic trap with frequencies 
uiy and u) z . The effective trap frequencies in y- and in¬ 
direction change by varying the potential of the lattice 
laser beams due to the change of the radial intensity 
profile. For the red-detuned lattice (red lattice) the fre¬ 
quency of the resulting additional confinement is related 
to its amplitude by 




/ \ V 2 
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and for the blue-detuned lattice (green lattice) the fre¬ 
quency of the resulting deconfinement is 


U a ,G(t) 


h (V G (t)\ 1/4 
mw a , G A G V E r ^c ) 


(4) 


Here a denotes the y- and ^-direction and w OJ is the 
waist of the corresponding lattice beam. The frequency 
in x-direction remains almost constant. 

The resulting time dependence of the total effective 
trap frequencies are 


WqW = [<4 i0 + W a,flO) - W a,G( 

^x {t) = ^x,0- 


, 1/2 


a = y,z 


(5) 


In order to be able to compare directly to the experiment 
[5(1 ] we use the initial trapping frequencies in x-, y- and z- 
direction u) x ,o = 2-7T x 38 Hz, u ) Vt o = 2-7T x 66 Hz and w 2j o = 
2-7T x 376 Hz. The wavelength of the infra-red laser beam 
is A r = 1064nm with waists of w 9j fl = 120/im, w Zi ^j = 
130/im and the wavelength of the green laser beam is 
A g = 532nm with the waists of w s G = 189/rm, w~ jG = 
42/rm in y- and z-direction. 
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FIG. 1. (Color online) The time dependence of the ampli¬ 
tudes of a) the red and b) the green lattice potential during 
the loading process. Insets: The black curves are the local po¬ 
tential for ~^ G < x < at a) t = 50ms when the red lattice 
is fully ramped up and b) t = 100ms when both the red and 
green lattice are ramped up. c) The evolution of the lattice 
potential V(x,t) shown for one double well during the load¬ 
ing process, d) The evolution of the local density distribution 
in one double well during the ramp up process. The density 
distribution is simulated for the homogeneous gas (case A). 


The single particle Hamiltonian (Eq. [T]) is separable 
in the different directions (H = H x + H y + H z ). The 
simulations are performed using exact diagonalization on 
the discretized version of the Hamiltonian and the full 
time evolution of the fermionic gas is determined. If not 
mentioned otherwise, the space is discretized using A* = 
^ in ^-direction, Ay = Az = in y- and z-direction 
and the time-step is discretized to At = 0.1ms. This 
value is chosen since it corresponds to the steps in which 
the experimental lattice is regulated in Ref. [56j . We 
discuss the convergence in these parameters later on. 


III. LOADING THE ATOMS INTO THE 
OPTICAL SUPERLATTICE 

In this section we investigate the loading procedure 
of the non-interacting fermionic gas into the optical su¬ 
perlattice potential along x-direction (Fig. [l]). We con¬ 
centrate hereby on loading the atoms first into the red, 
large wave-length optical lattice. The height of this lat¬ 
tice potential is changed slowly using an S'-shaped curve 
of the ramp. Only after the red lattice has reached its 
full height the short (green) wavelength optical lattice of 
wavelength A G is additionally ramped up using the same 
ramp form. The first part of the ramp therefore also rep¬ 
resents the typical loading procedure of Fermi gases into 
two-dimensional pancakes [5(| l58l - [62{ . 

The ramp sequence is shown in Fig. [lji.b. The optical 
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lattice potential in x-direction during the ramp process 
is described by Eq. [2] The intensity V 1 at time t follows 
the 5-shaped functional form 

V 7 (t) = V lti + (V 7l/ - V^i) sin 2 , (6) 

where 7 can be R and G for the red and green lattice, 
respectively. V~ /t i and V^j are the initial and final optical 
lattice potentials. The times t 7j o and r 7 are the start time 
of the ramp and the ramping time for the corresponding 
lattice and Vy(t) is valid for t 1 fi < t < t 7j0 + t 7 . The 
red lattice is ramped up in tr — 50ms from Vr^ = 0 to 

h 2 k 2 

Vrj = 48 E ri R where E t ,r = is the recoil energy of 
the red lattice (see Fig. [TJa.) . Afterwards the green lattice 
is ramped up for tq = 50ms from Vg,» = 0 to Vqj = 
11.7 E t ^g, where E rt o{= 4 E X} r) is the recoil energy for 
the green lattice (see Fig. QJd). In the loading process 
we fix the phase to <f> = 0.127T such that the final lattice 
potential consists of a strongly imbalanced superlattice. 
The evolution of the potential for a single double well is 
shown in Fig. [TJ3. 

In order to investigate the loading process we perform 
simulations starting with the ground state of a harmoni¬ 
cally trapped gas. We use the initial trapping frequencies 
in y- and z-direction as uJ y .o and uj z .o- In order to iden¬ 
tify the influence of the trapping potential along the x- 
direction we consider two different cases which we study 
in details in the two following sections: (A) a homoge¬ 
neous gas along the x-direction, i.e. u x = 0 and (B) a 
trapped gas with a weak trapping frequency co x ,o- 


A. Homogeneous gas 

In this section we consider the homogeneous setup 
with uj x = 0 and open boundary conditions along x- 
direction. Further, we choose the trap frequencies in y- 
and z-direction independent of time, i.e. u! y (t) = o>, y .o = 
27t x 66Hz and oj z {t) = ui Zj 0 = 27 r x 37 QHz. Here the 
space is discretized as Ax = ^. We calculate the den¬ 
sity distribution for 20 000 particles in 20 double wells (40 
single wells of the green lattice) during the ramping up. 
The evolution of the density distribution of a single dou¬ 
ble well is shown in Fig.Q]!. The initially homogeneously 
distributed density quickly (time « 10 ms) localizes at 
the bottom of the arising lattice wells of the red lattice 
potential. For the final red lattice depth of 48 E x r, the 
localization of the density is almost perfect and in be¬ 
tween the time 40ms and 50 ms only slight changes can 
be identified. When the short wave-length (green) lattice 
is ramped up starting at t = 50ms the localized density 
within one red lattice well splits into two parts. The 
dominating part of the occupation evolves into the deep 
odd lattice well of the superlattice, whereas only a very 
low occupation is transformed to the higher even lattice 
well. 

In order to monitor the excitations which are cre¬ 
ated during this ramping process, we determine the time 
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FIG. 2. (Color online) Time evolution of the instantaneous 
eigenenergies in x-, y- and z-direction during the ramping up 
process for the homogeneous case of 20 double wells with open 
boundary condition. The energy levels are shifted to have a 
minimum energy equal to zero at any time. The black lines 
are occupied levels and the red (grey) lines are unoccupied 
levels. 
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FIG. 3. (Color online) Time evolution of the occupation of the 
eigenenergy levels of H x integrated over the y- and z-direction 
for the case of the homogeneous system with 20 double wells 
and 20000 particles during the ramp up process. The vertical 
axis is the quasi-momentum index of the eigenenergy levels in 
x-direction and the first Brillouin zone of the red lattice lies 
between -1 and 1. 


evolution of the occupation of the instantaneous energy 
eigenlevels (the eigenenergies at each time t). Since the 
Hamiltonian does not change in time along the y- and z- 
direction and is separable, the eigenvalues of the parts H y 
and H z and their occupation remain constant over the en¬ 
tire ramp as can be seen in Fig.j^z. In contrast along the 
x-direction the ramp up of the lattice changes strongly 
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the energy levels (Fig. [2jc). During the turning on of the 
red lattice Bloch bands develop. Each band contains 20 
energy levels (due to the choice of finite length of 20 red 
lattice wells). The bands are well separated already after 
approximately 10ms. A drastic narrowing of the bands 
is observed. The switching on of the green lattice shifts 
the position of the bands but almost no mixing of the 
two Bloch bands is seen. In x-direction the evolution of 
the occupation of the energy levels depends strongly on 
the ramp procedure. For the chosen one, the evolution is 
shown in Fig. [2}c, where occupied levels are represented 
as black lines and in Fig. [3] The initial distribution due 
to Pauli blocking occupies around 26 lowest energy levels 
in x-direction and almost no change can be seen during 
the ramp. This means that the ramp procedure is almost 
adiabatic with respect to the x-direction. Comparing the 
final occupied energies in the different directions, one sees 
that even though the dynamics is almost adiabatic in re¬ 
direction, a highly excited final state is created. The 
ground state of the configuration at the end of the ramp 
would have lower occupation in x-direction and more lev¬ 
els occupied in y- and 2 -direction. This is caused by the 
decoupling of the different lattice directions. In an ex¬ 
periment a finite interaction or misalignment of lattices 
would cause a distribution of the fermionic occupations 
between the different lattice directions. Let us note that 
we need to be careful with the interpretation of our simu¬ 
lations at the initial times, since finite size and discretiza¬ 
tion effects might play an important role. We will study 
a large system subjected to a trap in the next section. 


B. Trapped gas 

In the presence of an additional trapping potential 
along the x-direction, the dynamics of the Fermi gas 
changes. We take the full time-dependence of the dif¬ 
ferent trap frequencies given by Eq. [5] into account. Here 
the time step used in the simulations is At = 0.5ms which 
is still much smaller than typical time scales of the in¬ 
duced changes beside at the initial time. Fig. 0] shows 
a zoom of the evolution of the density distribution dur¬ 
ing the ramping up. The density distribution starts from 
a Thomas-Fermi like profile at the initial time t = 0 
(cf. green curve in Fig. [5]). This continuous distribution 
is split by increasing the red lattice potential and after an 
initial ramp up {t < 20ms) the Fermions mainly occupy 
the minima of the red lattice wells. The actual density in 
a certain lattice well depends on the position of the lattice 
well in the trap (cf. t = 50ms line in Fig. 0] and the red V 
in Fig. [5j. In particular, the central wells contain more 
particles than the boundary wells, in which the density 
falls off to zero. The subsequent splitting of these red 
lattice wells into the superlattice configuration is almost 
adiabatic with respect to the local splitting. However, 
in the central wells, due to the initially higher density, 
the dynamics induces occupation in the even wells of the 
green lattice (cf. t = 100ms in Fig. 0] and the blue > in 
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FIG. 4. (Color online) Time evolution of the local density for 
N = 10 5 fermions with 250 double wells during the loading 
process (Case B). For better visibility only 5 double wells in 
the center of the trap are shown. 



FIG. 5. (Color online) The local density of the wells inte¬ 
grated over y- and 2 -direction and over the wells at different 
times of the ramping process for 10 5 fermions in the harmonic 
trap (case B) with 250 double wells. The ground state before 
and after the loading process are also shown. For each double 
well index there are two points showing the number of parti¬ 
cles in the even and odd well. The odd wells are deeper and 
thus, have higher density. The inset shows the local density 
inside one double well in the center of the trap. 


Fig. EJ). 

In the following we connect the dynamics to the change 
of the occupation of the instantaneous energy levels 
shown in Figs.[6]and0 As in the case of the homogeneous 
gas, the energy levels in x-direction respond strongly to 
the ramp up of the red lattice. Due to the time-dependent 
change of the trap frequencies in y- and 2 -direction, also 
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FIG. 6. (Color online) Evolution of the energy level occupa¬ 
tion in ^-direction integrated over y- and ^-direction during 
the ramping up process for N = 10 5 fermions in the har¬ 
monic trap (case B) with 250 double wells. The vertical axis 
is the quasi-momentum index of the eigenenergy levels in in¬ 
direction. Once the bands are separated, the first Brillouin 
zone of the red lattice lies between -1 and 1. 
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FIG. 7. (Color online) Time evolution of the instantaneous 
eigenenergies in x-, y- and ^-direction during the ramp up 
process for the harmonic trap (case B) with 250 double wells. 
The energy levels are shifted to have minimum energy equal to 
zero at any time. The black lines are occupied levels and the 
red (grey) lines are unoccupied levels for N = 10 5 fermions. 


the energy levels of H y and H, change. The change in 
y- and z-direction is found to be very smooth and the 
level occupation shows an adiabatic following of the lev¬ 
els in these directions. In contrast, along the ^-direction 
Bloch bands arise around approximately 10 ms, which 
separate clearly around 20ms. In contrast to the homo¬ 
geneous case, we see that the energy levels show a strong 
mixing which is due to the trap along x-direction. In 


FIG. 8. (Color online) Zoom of the time evolution of the 
instantaneous eigenenergies in x-direction during the ramping 
up process for the harmonic trap (case B) with 250 double 
wells. The energy levels are shifted to have minimum energy 
equal to zero at any time. Only 1 out of 10 energy levels is 
plotted. One can see the level crossings at the beginning of 
the ramp process. The inset shows the avoided level crossing 
of five neighbouring energy levels. 


particular, in Fig. [7] we see that occupied (black) levels 
arise from the lowest bands and split into Bloch bands 
around 10ms (Here we show three of the Bloch bands). 
If the ramping up was adiabatic, then we would expect 
to have all final states in the lower levels. However, there 
are many avoided level crossings (Fig. EJ in the begin¬ 
ning (before 10ms) and the system is excited to higher 
levels (Fig. Q. This can be seen well in the occupation 
of the instantaneous energy levels in Fig. [G] An adia¬ 
batic ramping would mean no change in the occupation. 
However, we see that a strong change of the distribution 
is seen up to approximately times of the order of 40ms. 
Thus, we see that in the trapped system the dynamics is 
strongly non-adiabatic caused by the non-locality of the 
eigenstates. The ramp up scale is too fast in order to 
assure a redistribution of the entire atomic cloud. 


IV. RELAXATION DYNAMICS DURING AND 
AFTER A QUENCH 

In this section we study the evolution of the prepared 
imbalanced gas after a sudden or slow ramp down of the 
red lattice. This change often called a quantum quench, 
sets the system out of equilibrium and induces a dynam¬ 
ics. 

We start with the initial patterned density distribution 
with alternating empty and occupied wells in the dichro¬ 
matic superlattice generated by the loading process de¬ 
scribed in section Mil By switching off the red lattice 
potential, the double well lattice structure turns into the 
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FIG. 9. (Color online) The time dependence of the red lattice 
potential during the slow quench (0 < t < 6ms). The black 
curves inside are the local potential for < x < -4f- before 
and after the quench. For sudden quench the ramp shape is 
the same but with smaller ramp time r = 0.5ms. 


normal green optical lattice. We assume the amplitude 
of the red lattice ramps down to zero with the same form 
as given in Eq. [ 6 ] (cf. Fig. 0. Two different ramp times 
will be considered (i) r r = 0.5ms for an almost abrupt 
quench which is faster than the tunneling times of the 
double wells and (ii) r r = 6 ms for a slow quench which is 
much larger than the final tunneling times in the green 
lattice. Some properties of the subsequent relaxation dy¬ 
namics following these quenches have been analysed by 
us and collaborators in Ref. [56j j . We complement these 
results here by a more detailed analysis showing addi¬ 
tional results. 

In Fig. [TO] the time evolution of the instantaneous 
eigenenergies in x -, y- and 2 -direction are shown. The 
black lines are the occupied levels during the slow quench. 
In x-direction the lowest two bands join at the end of the 
ramping down of the red lattice potential. The higher 
bands remain well separated from the two lower bands 
since they correspond to higher Bloch bands of the green 
lattice potential which is still present. During the last 
stage of the ramp down, many level crossings occur in 
between the lowest two bands. 

In Fig. [Til l and b the time evolution of energy level 
occupation is plotted for the almost abrupt and slow 
quench, respectively. In both cases the occupation re¬ 
mains almost constant over a large percentage of the 
quench time. Only during the last third of the quench 
time changes in the occupation are evident. This is due to 
the initially large gap in between the Bloch bands which 
prohibits tunneling. Only when the lowest two Bloch 
bands join, the occupation changes and a redistribution 
occurs. The final distribution of the occupation after the 
different quenches is shown in Fig. [12] In the case of 
the almost abrupt quench many excitations are created. 
This is reflected in the broad distribution created up to 
relatively large levels. In comparison during the slow 
quench, more fermions are transferred to the low energy 
levels. This process is signalled by the clear peak in the 
energy level occupation around zero. However, surpris¬ 
ingly the distribution shows a double structure and the 
number of excitations of the high energy levels (> 0.5) 
are comparable to the number of excitations in this re¬ 
gion generated by the almost sudden quench. 

In the following we study the relaxation dynamics of 
the system at different time scales after the almost abrupt 



FIG. 10. (Color online) Time evolution of the instantaneous 
eigenenergies in x-, y- and 2 -direction during the slow quench 
for the harmonic trap with 250 double wells. The energy levels 
are shifted to have minimum energy equal to zero at any time. 
The black lines are occupied levels and the red (grey) lines are 
unoccupied levels for N = 10 5 particles. 
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FIG. 11. (Color online) Time evolution of the energy level 
occupation in ^-direction integrated over y- and 2 -direction 
during (a) abrupt quench and (b) slow quench. N = 10 s 
fermions in the harmonic trap with 250 double wells. The 
scale of the vertical axis is chosen as in Fig. [6] 


and slow quench starting from these highly excited states. 


A. One-dimensional tight binding model 

Before presenting numerical results for the dynamics 
after the quench in three dimensions we start with a toy 
model for fermions in one dimension. The Hamiltonian 
is a tight-binding model with tunneling amplitude J and 
energy offset A between neighbouring sites 

L A 

H = 5>Jc}c, + i - Jc] +1 Cj + (-1 Y—c)c]} (7) 

l=i 

where L is the number of sites. This model has previously 
been investigated ffl ITtI . Mm for bosonic atoms. For 
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FIG. 12. (Color online) The energy level occupation in x- 
direction integrated over y- and ^-direction before and after 
the abrupt and slow quench. The scale of the horizontal axis 
is chosen as in Fig. [6] 

completeness we summarize some of the results in the fol¬ 
lowing. For A = 0, the energy band is the cosine-shaped 
band. At finite A a gap opens and splits the energy band. 
In the limit of y —> oo and half filling (N = -j) only odd 
sites are occupied by one fermion each. For a smaller 
energy offset, this imbalance decreases until for the ho¬ 
mogeneous lattice without energy offset, a homogeneous 
density distribution is reached. 

The system can be excited by a sudden quench of the 
energy offset from a finite value A to zero. The subse¬ 
quent time evolution of the local density rij = (CjCj) is 
determined by 

, N 1 , /4 Jt 2t rifcV , 

n j (t) = 2 - (- 1 y-j-J2 tajnh6k cos (— cos — 7 —) ( 8 ) 

^ k= 1 ' ' 

where sinh0fc = A 27rk . Thus, this dynamics depends 

4J|cos L | 

on the initial energy offset by the function tanh^.. 

We define the averaged odd-even contrast as C oe = 
N °~ N Ne where N 0 and N e are the total particle numbers 
in odd sites and even sites, respectively. This quan tity 
can be measured in the experiment 0, lU [ 50 L l56l l63j 
(section IV]) . In the limit of a large initial energy offset 
(y —> 00 ) and infinite size lattice, the averaged odd-even 
contrast is C oe (t) = Jo(^) where Jq is the zeroth-order 
Bessel function of the first kind (same as what found for 
bosons [Kl, El). The averaged contrast is maximum in 
the beginning C oe ( 0) = 1, then relaxes symmetrically 
about C oe = 0 with the period of T = yy and decays 
asymptotically in time as f 1 / 2 . If y is finite the initial 
odd-even contrast is less than 1 and it follows similarly 
damped oscillations than the Bessel function in time. It 
differs mainly in a lower amplitude and the details of the 
damping. 


Note that this tight-binding model only covers the po¬ 
tential imbalance and not the exact form of the double 
well which might include a dimerized hopping. Further, 
higher Bloch bands are completely neglected and a pure 
one-dimensional system without additional trapping po¬ 
tential is assumed. Nevertheless, we will recover some of 
the features in the full relaxation dynamics. 


B. Short time dynamics: Relaxation of the local 
imbalance 

In this section we determine the time evolution of the 
density distribution of the full model. By integrating 
the density distribution over y- and z-direction and also 
inside each well i we calculate the number of particles rq 
in the two dimensional sheet corresponding to the well i. 
Similar to the tight-binding model we probe the density 
imbalance of the odd-even wells An; = ri 2 i+i—n 2 i (where 
l is the double well index) and the averaged odd-even 
contrast C oe which is averaged over the trap. 

Fig. ITTil i shows the evolution of the odd-even density 
imbalance An; rescaled by no(t = 0) after the almost 
abrupt quench. The central well closely follows the ex¬ 
pected Bessel-like function of the one-dimensional tight- 
binding model with y = 4 in Eq. [H] (Fig. fl4h). The 
period in the Bessel function is T = where Jcj is the 
tunneling amplitude of the green lattice. We determine 
the tunneling amplitude from a band structure calcula¬ 
tion (Jg = 1.3085 x 10 ~ 2 E r) c)- The resulting period is 
T « 1.08ms which has slight deviation of 5% from the 
period obtained in the numerical results of the full model. 
This slight discrepancy is due to the chosen spatial dis¬ 
cretization of Ax = 4# and becomes smaller for smaller 
discretizations. The oscillations of the local odd-even 
imbalance at the boundaries of the cloud (at l « ±70) 
show a considerable slowing down and Bloch oscillations 
appear. The spreading of the period of the oscillations 
leads in the trap-averaged odd-even contrast to a stronger 
damping in time (the blue curve in Fig. ITTh). Both the 
slowing down of the tunneling and the onset of Bloch 
oscillations are due to the local energy offset between 
neighbouring wells caused by the trapping potential. The 
position of the maximum contrast closely follows the lo¬ 
cal Bloch period given by Tg(Z) = —tttti using the local 
offsets caused by the trap. The slight oscillating features 
along the x-direction are due to the initial non-adiabatic 
preparation of the atoms in the superlattice. As can be 
seen in Fig. 0 the density distribution after the quench 
has a low occupation in the upper Bloch band, which 
shows strong oscillatory structures along x-direction. 

In Fig. mb the time evolution of the odd-even den¬ 
sity imbalance An; rescaled by no(t = 0) after the slow 
quench is shown. The general features of the relaxation 
dynamics after the slow quench are very similar to the 
almost abrupt quench at short times. In particular, a 
damped oscillatory behaviour is observed. The main dif¬ 
ference is the initially lower amplitude. This can be bet- 











t (ms) t (ms) 


FIG. 13. (Color online) Evolution of the local odd-even den¬ 
sity imbalance rescaled by its maximum value after (a) the 
almost abrupt quench and (b) the slow quench for N = 10 s 
fermions in the harmonic trap with 250 double wells. The 
dashed lines show the predicted Bloch oscillations. 


ter seen in Fig. 114b . where the time evolution of trap- 
averaged odd-even contrast during and after the slow 
ramp down is plotted. The odd-even contrast at the be¬ 
ginning of the ramp (t = —6ms) is 0.92 which is the 
contrast after initial non-adiabatic preparation. For the 
almost abrupt quench the value at time t = 0ms is very 
close to this value. In contrast, during the slow quench, 
the odd-even contrast is much lower and the following os¬ 
cillations have a much smaller amplitude which is quickly 
damped out. The different behaviour is related to the 
more efficient transfer of the fermions in the last stage 
of the quench to low energy levels. These are mainly the 
energy levels which correspond to the local redistribution 
of the fermions, whereas higher energy levels which are 
still excited almost as intensively as in the abrupt quench 
correspond to extended states. Thus, the width of the 
density distribution has not yet varied considerably dur¬ 
ing the slow quench, such that no global relaxation has 
taken place. This motivates us to investigate in more de¬ 
tail the global relaxation taking place at longer times in 
the next section. 


C. Breathing modes 


Additionally to the dynamics along the x-direction, the 
change of the effective trapping frequency in y- and z- 
direction (Eqs. [5]) due to the increase of the lattice along 
x-direction induces oscillations of the density distribu¬ 
tion. In Fig. [T5] b) we show the local density along y- 
direction in the center of the cloud (x = z = 0) at dif¬ 
ferent times. We perform the simulation for a discrete 
tight-binding model as explained in the appendix [A] A 
breathing of the density distribution with a breathing pe¬ 
riod T « 10ms (corresponds to the breathing frequency 
2wj, « 27 t x 96 HZ) can be clearly seen. This strong 
change in the density distribution is also reflected in the 
evolution of the distributions along the other directions 
as seen in Fig. [15] a and c. Let us note, that this could 
not be seen in the previously discussed observables as the 


(a) 



FIG. 14. (Color online) a) Evolution of the local even-odd 
contrast (ci = ] in the center of the trap (l = 0), 

the boundary of the trap (l = 75), and the trap averaged 
contrast Coe after the almost abrupt quench. The contrast in 
the center fits well with the Bessel-like function predicted from 
the tight-binding calculation for the homogeneous system, b) 
Evolution of the trap averaged contrast C oe during and after 
the slow quench. For both cases the quench ends at t = 0. 



t (ms) 


FIG. 15. (Color online) Dynamics of the local density in y- 
direction (x = z = 0) for the discrete tight-binding model (see 
Appendix [A] after the abrupt quench for 10 5 particles in the 
trap with 500 sites. 


local density imbalance, since these were integrated over 
the y- and z-direction. 
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D. Long time: Global relaxation dynamics 

In this section we focus on the global relaxation dy¬ 
namics of the cloud after the quench. This is motivated 
by our finding in the previous section that the local odd- 
even density imbalance relaxes algebraically fast to zero, 
a possible thermal equilibrium value. The initial density 
profile along x —direction before the quench is relatively 
broad since only every second lattice well is occupied. 
The ground state profile in the lattice after the quench 
would be much more narrow. In the presence of the trap 
the global dynamics after the quench acts to redistribute 
the density profile. A good candidate to probe such a 
global relaxation inside the trap is the quasi-momentum 
distribution of the corresponding released system mea¬ 
sured in experiments by the band-mapping technique. 
In the simulations we obtain the quasi-momentum dis¬ 
tribution by a direct projection onto the homogeneous 
quasi-momenta. This distribution would not evolve with 
time in the absence of a trapping potential and is, thus, 
very well suited to measure the induced global relaxation. 

In Fig. M the time evolution of the homogeneous 
quasi-momentum distribution is shown after the slow 
quench. An initial focusing of the occupation towards 
the center (t < 10 ms) is observed for the quasi-particles 
in the first Brillouin zone ( \q\ < kc)- Subsequently, they 
are accelerated to higher momenta which is signalled by 
the outgoing diagonals from the center. The oscillations 
between low and high momenta which leads to the rep¬ 
etition of such diagonals, make the evolution resemble 
a fish-bone-like structure. A small fraction of particles 
are in the higher bands (« 10%). In the second Bril¬ 
louin zone (ko < |qj < 2 kc) a similar acceleration is seen 
which takes place with a different frequency. In both 
Brillouin zones, these oscillations seem to persist over a 
long time (several hundreds of ms). For the case of the 
almost abrupt quench the main feature resemble those 
of the slow quench. However, they are much less pro¬ 
nounced (Fig. [17]) and faster oscillating structures occur. 

In order to quantify the dynamics better we define the 
contrast i wit^rin the first aqd second Brillouin zones as 
Ci = and C 2 = where Nj is the num¬ 

ber of atoms in the first half of the j —th Brillouin zone 
i.e. the occupation of the quasi-momentum states in the 
interval < M < aR d A? is the number 

of atoms in the second half of the j—th Brillouin zone 
i.e. the occupation of the quasi-momentum states in the 
interval < M < ^ 7 . 

In Fig.fTBh the evolution of the contrast within the first 
Brillouin zone after the almost abrupt quench and slow 
quench is plotted. The contrast for the almost abrupt 
quench is scaled by a factor of « —2 and shifted by 
0.6ms. One sees that the qualitative behaviour of the 
contrast for both is very similar in particular the oscilla¬ 
tion periods are comparable (even though the contrast for 
the almost abrupt quench is much smaller). We observe 
damped oscillation with a period of T « 43ms. This pe- 



t (ms) 


FIG. 16. (Color online) Evolution of the homogeneous quasi- 
momentum distribution in i-direction integrated over y- and 
2 -direction after the slow quench. It is simulated for 10 5 par¬ 
ticles in the trap with 250 double wells. The scale of the 
vertical axis is chosen as in Fig. [S] 

riod is close to the breathing mode of the particles in the 
trap in the presence of the green lattice. The breathing 
period corresponds to the new rescaled trap frequency 

« ttlo x is T 2 ~ 36.6ms. 

In Fig. [T5b the time evolution of the contrast within 
the second Brillouin zone of the green lattice after the al¬ 
most abrupt quench and slow quench is plotted. The time 
is shifted by 2ms for the slow quench. The oscillating pe¬ 
riods for two different quenches are very close. Since the 
particles in the second Brillouin zone are living in the sec¬ 
ond Bloch band of the green lattice, the tunneling ampli¬ 
tude between neighbouring sites J G in the higher band is 
different. The breathing mode for particles in the second 
energy band corresponds to the new rescaled trap fre¬ 
quency IT ss . With the tunneling amplitude 

calculated for the second band ( J' G = 1.9492 x 10 _1 iH ri G) 
the breathing period is T 3 ps 9.5ms that is close to the 
oscillating period of the contrast in the second Brillouin 
zone. 


E. Discussion of the steady states 

In this section we investigate the steady states which 
arise at long times after the abrupt and slow quench. In 
Fig. [T 2 ] the energy level occupation after the quenches 
has already been shown. Since this distribution does not 
change anymore in the time-independent situation after 
the quench, this distribution corresponds to the steady 
state. We had seen that both, the distribution after the 
sudden and after the abrupt quench had excitations up to 
relative high energy levels. Further, the distribution after 
the slow quench has shown a distinct two peak structure 
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t (ms) 

FIG. 17. (Color online) Evolution of the homogeneous quasi¬ 
momentum distribution in i-direction integrated over y- and 
2 -direction after the almost abrupt quench. It is simulated 
for 10 5 particles in the trap with 250 double wells. The scale 
of the vertical axis is chosen as in Fig. [S] 

(a) 



FIG. 18. (Color online) Evolution of the quasi-momentum 
contrast within (a) the first Brillouin zone and (b) the second 
Brillouin zone after the abrupt and slow quench. For the 
almost abrupt quench in (a) the contrast (Ci) is scaled by a 
factor of approx. —2 and then shifted by approx. —0.6. For 
the slow quench in (b) the time is shifted by 2ms 


of a relatively high peak around low energies and a much 
broader peak with lower amplitude which extended up 
to high energies. In this section, we would like to see 
how this influences local observables as the density dis¬ 
tribution. The local density of the steady states after 
the slow and sudden quench are shown in Fig. [19] Addi¬ 
tionally, we display the density distribution after a com¬ 
pletely adiabatic quench. In both distributions, locally a 
balanced density has been reached which corresponds to 



FIG. 19. (Color online) Density distribution of the stationary 
state after the slow quench, the almost sudden quench, and 
the adiabatic quench. The density is integrated over y- and 
2 -direction and over the wells. 

the expected value in thermal equilibrium. However, the 
global density distribution shows similar features than 
the energy level occupation. After the slow quench, the 
global density distribution can be decomposed into two 
superimposed shapes. In the center the density distri¬ 
bution resembles the stationary state after an adiabatic 
quench (scaled by 0.73). Thus, in this region the density 
seems to come (at least from its form) close to a thermal 
equilibrium distribution. However, at the boundaries the 
density distribution resembles the stationary state after 
the sudden quench (scaled by 0.89) and its form is very 
distinct from a thermal distribution. We verified this 
by comparing the density distribution after the abrupt 
quench to a thermal distribution varying the tempera¬ 
ture. Therefore, the global distribution of the local den¬ 
sity signals strongly the non-equilibrium nature of the 
state. 


V. EXPERIMENTAL DETECTION 
PROCEDURE OF THE ODD-EVEN DENSITY 
CONTRAST 

In order to measure experimentally the initial relax¬ 
ation discussed in section IIV Bl the even-odd contrast 
has been monitored. The experimental detection of the 
even-odd contrast relies on the use of the bichromatic 
lattice IlHii & The experimental sequence 
consists of three parts (see Fig. [20]) : i) The red lattice is 
quickly {t\ = 0.1ms) ramped up to Vr = 48 E t ^r (at a 
phase (f> = -|) which produces a superlattice with strong 
energy offset between neighbouring sites. By this, the 
occupation of the green lattice well is projected onto 
the occupation of the corresponding well of the super¬ 
lattice. ii) The green lattice is ramped down to zero in 
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FIG. 20. (Color online) Ramp schedule of the depth of a) the 
red and b) the green lattice potential during the detection 
procedure of the integrated odd-even contrast. Time zero is 
chosen to correspond to the beginning of the ramp down of 
the green lattice. 



FIG. 21. (Color online) a) Time evolution of different frac¬ 
tions of the atoms in different Brillouin zones. The shaded 
area marks the time regime when some atoms reach the 
boundary and the simulation does not correspond to ex¬ 
perimental results any longer, b) Evolution of the quasi¬ 
momentum distribution in ^-direction integrated over y- and 
^-direction for early times of the measurement procedure. The 
scale of the vertical axis is chosen as in Fig. [ 6 ] c) The con¬ 
trast Ci/ 3,4 after the green lattice is ramped down during the 
measurement for different durations of the time steps. The 
simulations are performed for N = 10 5 particles in a trap 
where 450 double wells are taken into account. 


t ~2 = 0.5 ms. The parameters for the superlattice are cho¬ 
sen, such that during this ramp down the population of 
the even wells transfers to the third and fourth Bloch 
band of the red lattice and the atoms in the odd wells 
remain in the lowest band. The second band remains 
empty, iii) Subsequently an adiabatic band mapping is 
performed in order to measure the percentage of atoms 
in the different Bloch bands. This is performed by a slow 


ramp down of the red lattice in 73 = 6 ms with an ad¬ 
ditional time-of-flight imaging of the atoms. Thus, the 
measured time-of-flight images reveal the odd and even 
lattice well population integrated over the trap. The odd- 
even contrast ( C oe ) can be measured by the new contrast 
defined as C\j 3)4 = w ^ ere ^3 the number 

of atoms in the j -th Brillouin zone. 

In order to analyze imperfections of the experimen¬ 
tal detection process, we perform numerical simulations 
starting from the ground state of the green lattice (Vg = 
11.7 E r>G , V R = 0 , 0 = j), i.e. a balanced density config¬ 
uration ( C oe = 0). The quasi-momentum distribution is 
calculated at each time step by a direct projection onto 
the homogeneous quasi-momenta. 

In Fig. l2ll i the results of the theoretical simulations 
of the fraction of atoms in different Brillouin zones and 
the contrast Ci/ 3,4 during the measurement procedure 
are shown. The shaded area marks when the particles 
reach the boundary and our simulation is not reliable 
anymore. For the simulation, a time step of At = 0.1 ms 
is chosen, since this corresponds to the regulation steps 
of the experiment [Hj- Surprisingly, this discretization 
time-step causes an increase of the contrast during the 
ramp down of the green lattice. After the ramp down 
the contrast remains almost constant. The sharp increase 
is mainly due to the excitation of atoms from the third 
and fourth bands to the fifth and sixth bands as can be 
clearly seen in Fig. l2Tb . This increase of the contrast 
falsifies the detection outcome. In the considered case 
a contrast of Ci/ 3,4 ~ 0.11 would be measured instead 
of the expected zero value of the prepared state. We 
repeated the simulations for different initial contrasts. 
The effect of the measurement procedure is strongest on 
the shown initially balanced situation. On the totally 
imbalanced situation a discrepancy of about 3% is found. 
The slight shift of the zero contrast to a final value has 
been observed in Ref. [5(|. This shift could be avoided 
by decreasing the regulation step of the ramp down of 
the lattice potential. In order to verify this, we have 
performed simulations with varying time-step between 
At = 0.1 ms to 0.01 ms and we find that the fraction of 
atoms excited to the higher bands decreases to Ci/ 3,4 ~ 
0.01 (see Fig. l2Tb). 


VI. CONCLUSION 

In this work we have studied theoretically the load¬ 
ing procedure of non-interacting fermions into a one- 
dimensional super lattice, the relaxation after an abrupt 
change of the superlattice, and the measurement proce¬ 
dure of the odd-even density contrast. Our calculations 
have taken into account the full three dimensional and 
spatially inhomogeneous nature of the cloud subjected to 
an optical superlattice potential at the expense of con¬ 
sidering non-interacting fermions. The findings we have 
obtained led to insights which can be used to improve 
the experimental loading and detection procedures. Con- 
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cerning the loading ramp of the fermions into the optical 
lattice, in our results a large mixing and the creation 
of excitations has been observed only in the very initial 
time of the loading procedure until the Bloch bands be¬ 
come separated (even in the presence of a trap). Thus, 
the initial loading has to be optimized in an experiment, 
whereas after the opening of the gap a faster ramp can be 
applied without harm. Additionally for fermionic atoms 
-in contrast to bosonic atoms- the number of atoms and 
the trap configurations plays an important role, since 
due to the Pauli blocking each available energy level can 
only be occupied by one fermion (or two for the spinful 
situation). Starting from a certain density the higher 
bands have to contribute and a loading of these can¬ 
not be avoided even by a locally adiabatic loading pro¬ 
cedure. The solution to this problem is either to start 
with a low density or enough available levels in the di¬ 
rections perpendicular to the lattice direction. The trap 
has the additional problem that due to the lattice ramp 
up a change of the global density distribution is induced. 
Another way of seeing the strong influence of the trap 
is the induced level crossing during the ramp. Control¬ 
ling the frequency of the trap this effect can be reduced. 
The redistribution of the density takes typically a long 
time. The redistribution at the boundaries can be almost 
blocked by Bloch oscillations due to the local energy off¬ 
set. In the considered case of non-interacting fermions 
this redistribution is even at infinite time prohibited by 
the conservation laws and it will never relax to its ther¬ 
mal equilibrium value. Even though, this is an artefact of 
the integrable situation, some of these integrable features 
might also dominate the situation in which the integra- 
bility is slightly broken by a weak interaction [64}. Let 
us mention that a very slow redistribution had already 
been identified previously, for the interacting case. In 
the interacting model Mott-insulating barriers can arise 
which can almost totally block the redistribution of the 
atoms [til . [66j . Thus in experiments this is a crucial 
point which has to be considered. A corresponding ad¬ 
justment of the trapping potential, as it is partially the 
case in red-detuned lattices [67}, could cure this problem 
(cf. also recent quasi-exact simulations for bosonic atoms 
in one-dimension which support this suggestion @)- 

The relaxation dynamics after a sudden and slow 
change of the bichromatic potential to a monochromatic 
potential has shown that locally a fast relaxation can oc¬ 
cur for some observables, whereas the global relaxation 
takes place on a much longer time scale. Due to the 
integrability of our system, the long-time steady state 
could show a structure of two superimposed peaks with 
different width, which clearly does not correspond to a 
thermal equilibrium state. We identified different impor¬ 
tant time-scales. The local relaxation has taken place 
on the time-scale of the tunneling time, here given by 
Tl = ~ l-08ms. Additionally, the global density re¬ 

distribution occurs on a longer time scale. An impression 
of these give the breathing modes, which are approxi¬ 
mately one to two orders of magnitudes larger. 


In the detection procedure for the odd-even density 
contrast also the time-scales play a very important role. 
We have found that a fast switch off of the green lattice 
potential and dominantly a relatively long regulation step 
of the lattice potential influences the resulting contrast. 
Whereas this change is negligible for a high contrast, in 
the situation of vanishing contrast the results are falsified 
and a finite contrast is found. This can be cured in the 
experimental results by a change of the lattice ramp time 
and the regulation step. 


Appendix A: Comparison of continuous simulations 
versus discrete simulations 


In this appendix we describe the simulations performed 
in the discrete tight-binding approximation shown in 
Fig. [15] For the strongly imbalanced superlattice along 
ir-direction, the Hamiltonian of the system is 


L/2 

H x = ^ [— Jl2C2l-l C 2l — J23ch c 2l+1 

l=-L/2 

— J 13^21 — \ C2l-\-1 ~ J2Ac\ l C2l+2 + h.C.] 


L 

£ 

j=-L L 


i-iycUj + VtfcUj 


(Al) 


Here we label each double well with the index l and the 
index for the wells inside the double well l are j = 21 — 1 
and j = 21. We restricted the description to include hop¬ 
ping terms up to the second nearest neighbours. Here 
J mn is the hopping amplitude from the well j = m to the 
well j = n. Vt is the trap potential in the discrete form 
and A is the energy offset between neighbouring sites in 
one double well without the trap. With the parameters 
we used for the continuous simulation the energy offset in 
the initially prepared state is given by A = 15 E r n which 
is much larger than the tunneling amplitudes. In the sim¬ 
ulations we calculate the ground state of the Hamiltonian 
where H x is as Eq. IA1I and H y and H z as in Eq. [T] with 
the corresponding trap frequency when both of the red 
and green lattice are on. 

After the quench when only the green lattice is on and 
there is no double well structure any more, the Hamilto¬ 
nian in ^-direction is given by 

L 

H x = 'y ' [— JcCjCj -(-1 — JGCj+i c j + Vf.j CjCj]- 

J=-L 

Here Jq = 1-3085 x 10 ~ 2 E r ^G is the hopping in the re¬ 
maining green lattice. The trap potential Vt is the same 
as before the quench. H y and H z are the same as given 
in Eq. [T] with the corresponding trap frequency for the 
green lattice. 

In Fig. [2H we compare the local density evolution of 
the continuous (Eq. |T]) and tight-binding model (Eq. lAll) 
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starting from the ground state of the imbalanced dichro¬ 
matic lattice before the quench. The evolution of the lo¬ 
cal density is plotted for different wells in the trap. The 
density is summed over the y- and ^-direction and for 
the continuous case the density is also summed over each 
well. One can see that the behaviour for the discrete 
and the continuous model are very similar only slight 
differences in the frequency are seen which start to cause 
deviations after a time of 5ms. This we attribute to ei¬ 
ther the discretization of the continuous model or the 
approximate hoppings (J mn ) calculated from the band 
structure. 


0 5 10 15 20 25 

t (ms) 

FIG. 22. (Color online) Evolution of the local density of odd 
wells at different regions of the trap after the abrupt quench 
compared for the tight-binding and continuous model. It is 
simulated for 10 5 particles in the trap with 250 double wells. 
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